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>  A  new  stochastic  dynamic  target  model  is  proposed  on  the  assumption  that 
certain  targets  execute  evasive  maneuvers  orthogonal  to  their  velocity  vector. 
Along  with  this  new  acceleration  dynamic  model,  the  orthogonality  is  also 
enforced  by  the  addition  of  a  fictitious  auxiliary  measurement.  The  target 
states  are  estimated  by  the  modified  gain  extended  Kalman  filter,  and  the 
angular  target  maneuver  rate  is  constructed  on-line.  A  guidance  law  that 
minimizes  a  quadratic  performance  index  subject  to  the  assumed  stochastic 
engagement  dynamics  that  includes  state  dependent  noise  is  derived.  This 
guidance  law  is  determined  in  closed  form  where  the  gains  are  an  explicit 
function  of  the  estimated  target  maneuver  rate  as  well  as  time  to  go.  The 
numerical  simulation  for  the  two-dimensional  angle-only  measurement  case 
indicates  that  the  proposed  target  model  leads  to  significant  improvement 
in  the  estimation  of  the  target  states.  Furthermore,  the  effect  on 
terminal  miss  distance  using  this  new  guidance  scheme  is  given  and  compared 
to  the  Gauss-Markov  model.  ; 
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Abstract 

A  new  stochastic  dynamic  target  model  is  proposed 
on  the  assumption  that  certain  targets  execute  evasive 
maneuvers  orthogonal  to  their  velocity  vector.  Along 
with  this  new  acceleration  dynamic  model,  the  orthog¬ 
onality  is  also  enforced  by  the  addition  of  a  fictitious 
auxiliary  measurement.  The  target  states  are  estimated 
by  the  modified  gain  extended  Kalman  filter,  and  the 
angular  target  maneuver  rate  is  constructed  on-line.  A 
guidance  law  that  minimizes  a  quadratic  performance 
index  subject  to  the  assumed  stochastic  engagement  dy¬ 
namics  that  includes  state  dependent  noise  is  derived. 
This  guidance  law  is  determined  in  closed  form  where 
the  gains  are  an  explicit  function  of  the  estimated  target 
maneuver  rate  as  well  as  time  to  go.  The  numerical  sim¬ 
ulation  for  the  two-dimensional  angle-only  measurement 
case  indicates  that  the  proposed  target  model  leads  to 
significant  improvement  in  the  estimation  of  the  target 
states.  Furthermore,  the  effect  on  terminal  miss  distance 
using  this  new  guidance  scheme  is  given  and  compared 
to  the  Gauss-Markov  model. 

1.  Introduction 

The  target  tracking  problem  for  homing  missile  guid¬ 
ance  involves  the  problem  of  estimating  large  and  rapidly 
changing  target  accelerations.  The  time  history  of  target 
motion  is  inherently  a  jump  process  where  the  acceler¬ 
ation  levels  and  switching  times  are  unknown  a  priori. 
Due  to  this  arbitrary  and  unpredictable  nature  of  tar¬ 
get  maneuverability,  target  acceleration  cannot  easily  be 
modeled. 

A  considerable  number  of  tracking  methods  for  ma¬ 
neuvering  targets  have  been  proposed  and  developed 
for  both  new  target  models  and  filtering  techniques[l]- 
[8].  In  spite  of  the  numerous  modeling  and  filtering 
techniques  available,  target  acceleration  estimation  us- 
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ing  angle-only  measurements  is  relatively  poor.  Usually, 
the  target  tracking  problem  is  approached  by  model¬ 
ing  target  acceleration  with  a  first-order  Gauss-Markov 
model  and  applying  the  extended  Kalman  filter(EKF). 
One  difficulty  with  the  Gauss-Markov  model  is  that 
the  assumed  large  process  noise  spectral  density  induces 
Kalman  filter  divergence  even  when  the  target  maneu¬ 
ver  is  not  present  and  can  lead  to  a  target  acceleration 
magnitude  estimate  which  exceeds  the  actual  maximum. 
Another  problem  is  that  the  target  motion  is  not  well 
represented  by  Gauss-Markov  diffusion  process. 

In  an  effort  to  alleviate  these  problems,  the  circu¬ 
lar  target  model  has  been  proposed  as  a  target  motion 
model,  where  the  phase  angle  is  a  Brownian  motion  pro¬ 
cess  and  the  acceleration  magnitude  can  be  either  a  ran¬ 
dom  variable  or  a  bounded  stochastic  process.  This  tar¬ 
get  model  was  suggested  in  [4]  using  concepts  extracted 
from  [9,10]. 

By  including  a  priori  knowledge  of  the  target  motion, 
improved  estimates  of  the  target  states  can  be  obtained. 
This  idea  is  included  in  [6,7]  by  using  a  target  accelera¬ 
tion  model  which  employs  a  target  mean  jerk  term.  For 
conventional  targets  such  as  winged  aircraft,  the  longi¬ 
tudinal  acceleration  component  is  often  negligible  com¬ 
pared  to  the  lateral  component  in  evasive  maneuvers. 
This  notion  fits  the  circular  target  model  where  the  an¬ 
gular  rate  term  is  estimated  to  account  for  the  actual 
dynamics  of  the  coordinated  turn.  This  model  is  pre¬ 
sented  in  Section  2.  However,  an  approximate  state  ex¬ 
pansion  is  required  to  handle  the  unknown  angular  rate 
in  the  target  model.  This  approximate  dynamical  sys¬ 
tem  used  for  estimation  is  presented  in  Sections  3.1  and 
3.2.  Furthermore,  in  Section  3.3  the  orthogonality  be¬ 
tween  target  velocity  and  acceleration  can  be  viewed  as 
a  kinematic  constraint  where  compliance  is  enforced  by 
including  this  constraint  as  a  pseudo-measurement[7,8j. 
The  approximate  target  dynamics  and  pseudomeasure- 
ment  are  included  in  the  modified  gain  extended  Kalman 
filter(MGEKF)[ll]  and  is  presented  in  Section  4.  The 
MGEKF  is  selected  because  of  its  superior  performance 
over  the  EKF  especially  for  bearing-only  problems.  In 
Section  5,  a  linear  quadratic  guidance  law  is  derived  for 
this  circular  target  model.  This  guidance  law  remains  a 
linear  function  of  the  estimated  states,  but  the  guidance 
gains  obtained  in  closed  form  are  a  nonlinear  function 
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of  the  estimated  rotation  rate  and  time  to  go.  Finally,  a 
numerical  simulation  is  performed  for  a  two-dimensional 
homing  missile  intercept  problem.  Both  the  estimation 
process  and  the  terminal  miss  are  enhanced  by  the  new 
models  and  the  associated  estimator  and  guidance  law 
in  comparison  with  the  Gauss-Markov  model 

2.  Target  Acceleration  Model 

In  this  section  the  circular  target  acceleration  model 
is  presented,  and  the  dynamic  consistency  between  this 
target  model  and  an  assumed  nonlinear  target  model  is 
discussed. 

2.1.  Circular  Target  Motion 

The  two-dimensional  homing  missile  guidance  sce¬ 
nario  is  described  by  two  sets  of  nonlinear  dynamic  equa¬ 
tions  of  motion  for  the  missile  and  target 
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where  (xm  ,  xx)  and  ( yM  ,  yx)  are  inertial  coordinates, 
Vm  and  Vx  are  the  velocities,  a\ft,  a\f„  and  ax  are  the 
accelerations,  and  6m  and  6x  are  the  flight  path  angles 
[Fig.  2].  The  subscripts  M  and  T  denote  the  missile 
and  target,  respectively,  and  aM,  and  a  at,  are  tangential 
and  normal  accelerations,  respectively.  Only  the  normal 
component  of  the  acceleration  contributes  to  changing 
angular  orientations  of  each  vehicle,  and  the  target  is 
assumed  to  fly  at  constant  speed. 

The  following  target  model  is  assumed  to  be  used  in 
the  filter.  The  objective  is  to  choose  a  model  that  is 
linear  in  order  to  reduce  the  numerical  computation  of 
the  filter,  but  reasonably  consistent  with  the  nonlinear 
model  so  that  the  estimates  are  of  good  quality.  The 
target  model  for  the  filter  in  two  dimensions  is 

aT*(0  =  aTcos(u>t  +  6),  ax>(<)  =  axsin(ut  +  6)  (2) 

where  ax  is  a  constant  which  is  unknown  a  priori,  u>  is 
the  angular  velocity  to  be  estimated  in  a  right-handed 
coordinate  system,  and  6  is  a  Brownian  motion  pro¬ 
cess  to  represent  random  target  maenuver  phasing  with 
statistics 

E[dO]  =  0,  E[d6 2]  =  Qdt,  0  =  l/r0  (3) 

Here,  Q  is  the  power  spectral  density  of  the  process  and 
r©  is  the  coherence  time,  the  time  for  the  standard  devi¬ 
ation  of  9  to  reach  one  radian.  While  in  the  previous  cir¬ 
cular  target  model[4]  the  acceleration  components  were 
just  a  diffusion  procesc  along  a  circle,  those  in  the  new 
model  are  related  to  the  actual  target  motion  through  a 
term  of  physical  meaning,  w. 

2.2.  Dynamic  Consistency 


In  order  to  see  how  the  current  model  approximates 
the  assumed  nonlinear  target  dynamics(l),  consider  a 
deterministically  equivalent  case.  Integration  of  (2)  with 
0  =  0  and  M  >  0  yields 


Q>T  (XT' 

Vx  =  — si  nut,  Vv  — - cosui 

w  u> 


VT+—.  (4) 

Uf 


where  the  initial  conditions  are  Vx  =  0  and  Vy  =  —  Vj. 
Adding  the  square  of  each  component  and  moving  all 
terms  to  the  ieft  hand  side  gives 

2(1 -cosw<)[(—  )2-—  V'T]  =  0.  (5) 

id  id 

For  this  equation  to  hold  for  all  t  >  0 


which  is  equivalent  to  the  differential  equation  for  angu¬ 
lar  rate  in  (1). 

Furthermore,  by  taking  the  dot  product  of  the  velocity 
vector  with  the  acceleration  vector,  we  obtain 

Vx.dx  —  ax[—Vx  H — — Jsinwt  =  0  (7) 

id 

using  (6)  for  all  t  >  0.  This  orthogonality  between  tar¬ 
get  velocity  and  acceleration  demonstrates  the  dynamic 
consistency  of  the  proposed  target  acceleration  model 
for  the  filter  with  the  nonlinear  target  dynamics. 

3.  New  Dynamic  and  Measurement  Models 
for  Estimation 

The  previous  section  dealt  with  a  new  circular  model 
for  filter  implementation  in  order  to  exploit  an  assumed 
characterization  of  the  motion  of  a  typical  target.  In 
this  section,  the  stochastic  dynamic  equations  for  the 
new  target  r  je,  c.re  derived.  Furthermore,  the  kine¬ 
matic  fictitious  -  .surement  suggested  in  [7,8]  is  also 
discussed. 

3.1.  Formulation  and  Approximation 

Ito  stochastic  calculus[12]  applied  to  the  Eq.(2)  results 
in  a  stochastic  differential  equation  with  white  state- 
dependent  noise[9]. 


where  the  elements  —  j  in  the  drift  coefficient  are  the 
Ito  correction  ternis.  Note  that  the  problem  is  nonlinear 
due  to  the  unknown  w.  To  avoid  solving  the  nonlinear 
problem,  u  is  approximately  removed  by  an  expansion 
of  state  variables  as  in  [10].  Define  r.cw  states  as 


a*  =  u>ax 
a\  =  u>ay 


lfl-1 


with  the  assumption 


M  =  I^H«i-  (io) 

By  augmenting  the  dynamics  of  these  new  states  to  (2), 
an  approximate  dynamical  model  which  includes  this 
new  target  model  is  truncated  as 


Note  that  u  does  not  appear  explicitly  in  (11),  and  (11) 
is  a  linear  stochastic  differential  equation.  An  idea  of 
this  sort,  gi%’en  in  [10]  for  a  scalar  problem,  led  to  sig¬ 
nificantly  improved  filter  performance. 

3.2.  Two-dimensional  Intercept 
The  two  dimensional  intercept  problem  is  developed 
in  a  relative  inertial  coordinate  system.  The  system  dy¬ 
namics  are  expressed  in  the  following  set  of  equations; 

=  ur,  yr  =  vr, 

ur  =  aTl_  -  aMxr ,  i)r  =  aTyr  -  aMyr , 

aT*r  —  arcos(uit  +  6),  ary  =  arsin(u/t  +  0). 

(12) 

W  ith  the  assumption  of  u  <  1  and  truncation  of  the 
target  dynamics  up  to  the  second  order,  the  ten-element 
state  vector  is  defined  as 

1  =  [*r  Vr  Ur  vr  ax  ay  a*  aj  a\  ay]T 

=  [Zi  X2  X3  *4  *5  *6  *7  Z*  *9  r10]r. 

mu  ■  .  (13) 

Thus,  in  terms  of  the  expanded  state  space,  the  linear, 

stochastic  state  differential  equation  is  described  by 


In  Ref.[7,8]  the  filter  performance  is  improved  by 
introducing  a  kinematic  constraint  based  on  a  priori 
knowledge,  which  is  implemented  in  the  form  of  an  aug¬ 
mented  fictitious  measurement.  In  particular,  the  accel¬ 
eration  vector  is  assumed  to  be  related  to  the  velocity 
as 

Vt  -  dj  —  0  (18) 

under  the  assumption  that  the  target  accelerates  pre¬ 
dominantly  orthogonally  to  its  velocity  vector.  When 
this  condition  is  not  met,  the  acceleration  has  a  compo¬ 
nent  in  the  direction  of  the  target  velocity  as 

Vt  St  =  i)  (19) 

where  Vt  and  St  are  assumed  to  be  random  vectors  rep¬ 
resenting  target  velocity  and  acceleration,  and  rj  is  the 
uncertainty  in  the  orthogonality.  This  idea  can  be  im¬ 
plemented  in  the  form  of  a  discrete  pseudo-measurement 
as 

z2(tk)  =  h2(x(tk))  +  TJk 

=  VT,raT,r  +  VTyr  aTyr  +  T)k  (20) 

where  %  is  a  white  random  sequence  with  assumed 
statistics 

E[r)k]  =  0,  E[r,ktf]  =  V26k,.  (21) 

Note  that  the  variance  of  the  measurement  noise  corre¬ 
sponds  to  the  tightness  of  the  constraint.  In  other  words, 
the  larger  the  variance,  the  more  relaxed  are  the  require¬ 
ments  on  longitudinal  acceleration.  This  fictitious  mea¬ 
surement  is  used  along  with  the  angle  measurement  in 
the  modified  gain  extended  Kalman  filter  described  in 
the  next  section. 

4.  Estimation  of  Target  States 


dx  -  (Fx  +  Bu)di  +  Gd6  (14) 

where  F  is  given  by  the  relations  in  (12)  and  by  the 
coefficient  matrix  of  x  in  (11)  where  i  =  2,  B  is  a 
10  x  2  matrix  of  zero  except  for  £3!  =  S42  =  —  1, 
u  =  ( aM,r,aMVr)T ,  and 


G  =  [0  0  0  0  -  i(  rs  -i8  i7  -  no  xs  ]T.  (15) 

•3.3.  Measurement 
Angle  measurement 

Angle  information  in  discrete-time  is  assumed.  Then, 
the  measurement  at  time  tk  is 


F(tk)  =  hi{x{tk})+vk  =  tan~\yr/xr)+vk  =  x*(t*)+r* 

(16) 

where  vk  is  a  white  random  sequence  with  statistics 

E{vk}  =  0,  E[vkvJ]  =  \\6ki-  (17) 


In  this  section,  the  modified  gain  extended  Kalman 
filter(MGEKF)[ll]  is  derived  for  the  circular  target 
model  and  for  the  fictitious  and  angle  measurement  de¬ 
fined  in  the  previous  section.  Also,  considered  is  the 
method  to  reconstruct  the  maneuver  rate  using  the  esti¬ 
mated  states.  Given  the  continuous-time  dynamics  and 
discrete-time  measurements,  as  in  the  previous  section, 
construction  of  the  filter  is  completed  by  specifying  the 
time  propagation  and  measurement  update  procedures. 

4.1  Time  Propagation 

The  state  estimate  x(< /<» _  1  )is  propagated  from  the 
current  time  t,_j  to  the  next  sample  time  t,  by  integrat¬ 
ing 

HU1-:  :}  ~  F*(t/ti-i)+  Ru(t). 

h/U-i)  =  FP(t/U-1)  +  P(t/ti^l)FT  +  E[GQC?], 

X{t)  =  FX{l)  +  X{t)FT  +  E{GQGT). 
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(22) 


given  the  posteriori  estimate  x(<,_i/<j_i)  =  x(U_i)  and 
posteriori  pseudo-error  variance  P(t<_i/tj_i)  =  i). 

The  notation  /?(</<;_  i)  denotes  the  value  of  some  quan¬ 
tity  R  at  time  i  given  the  measurement  sequence  up  to 
time  <,■_  j.  The  integration  of  the  covariance  of  the  state 
X(t)  begins  with  A'(0)  at  time  t  =  0.  Upon  integrat¬ 
ing  the  equations  above  to  the  next  sample  time,  the 
propagated  estimates  are  obtained  as  follows: 

x(U)  =  x(U/ti- ,),  P(ti)  =  P(U/U- 1)  (23) 


It  should  be  noted  that  since  the  process  noise  is  state- 
dependent  ,  the  integration  to  propagate  P(t)  also  re¬ 
quires  the  integration  of  A'(/),  where  E[GQGT]  matrix 
turns  out  to  have  nonzero  elements  for  its  lower-right 
6x6  matrix.  Note  that  the  approximation  technique 
reduces  the  originally  nonlinear  dynamics  to  linear  dy¬ 
namics.  This  allows  for  the  closed  form  solutions  of  the 
propagation  of  the  estimates  rather  than  performing  on¬ 
line  integration. 

4.2  Measurement  Update[ll] 

The  states  are  updated  as  follows: 


A'(fi) 


where 


=  x(ti)  +  K(ti)[z-h(x(U))], 

=  P(ti)H(ti)T[H(ti)P(ti)HT(ti)  +  V]~l 

(24) 


V  = 


Vi  0 
0  v2 


and 


n>.) 


m) 


Hn, 

0, 


o.  o,  ,0,  ... 

Hz>,  Hlt,  H2i,  Hm, 


(25) 

:] 


with 


Hu 

Hu 


~T~^,  Hu 

+  y, 

Is,  Hu  =  it. 


i?  +  y?  ’ 

Hu  —  X3  +  >  HK  =  x«  + 


where  the  missile  acceleration  in  the  x  and  y  directions, 
a\it  and  are  assumed  to  be  measured  very  accu¬ 
rately  with  on-board  sensors. 

The  measurement  update  of  the  pseudo-error  variance 
is  performed  by 


PM 


[/  -  K{U)g{z(ti),x(U 


(26) 


T 


where  g(z(ti),x(ti))  is  used  in  the  update  of  P  rather 
than  H  of  Eq.  (25)  and  is  given  as 


h(x{t,))  -  h(x(ti))  = 


tan 


Ir(ij)  Xr(ti) 


=  g(z(ti),£(U))(x(U)  -  *(<.))  (27) 

Note  that  g  is  a  2  x  10  matrix  of  function  explicit  only 
in  the  known  quantities  z  and  x.  In  this  sense,  the  func¬ 
tion  h  has  a  universal  linearization  with  respect  to  the 


measurement  function  z.  Unfortunately,  this  type  of 
linearization  with  respect  to  the  measurements  occurs 
for  only  a  few  functions.  It  is  applicable  to  angle  mea- 
surements[ll],  but  not  for  our  new  pseudo-measurement. 
Therefore,  we  must  for  the  pseudo-measurement  revert 
back  to  the  extended  Kalman  filter  form  and  define 


9i{*(U),x(ii))(x(U)  ~  x(U))  =  Mx(<<))  “ 

dh2 


dx(ti) 


|r=*(I(t|)  —  x(ti)) 
(28) 


dh7 


where  the  expression  for  „  .  .  . 

dx{U) 

ond  row  of  H  in  (25) 

For  angle  measurements^  1] 


-t  is  found  in  the  sec- 


M*(*0)  -  Ai(*(*i))  =  0i(*(<<).  *(<<))(*(*<)  -  *(*<)) 

=  -  E(ti)H(z*(ti))(x(t;)  -  x(U)) 

(29) 


where 


E(U)  = 
D(t$)  = 


D(ti)tan~1a(ti) 

a(tj) 

yJxT(tt)'2  +  yr(tj)2 


( Xr(ti)Xr(ii )  +  3/r(t;)yr(t|)) 

a(U)  =  Z>(<,)tf(z*  (!,))*(*< ) 

H(z*(ti))  =  [sin  z*(U)t  -cos  z'(tj)  ,0,0,0  ,...] 


(30) 


As  discussed  in  [11],  g(z(t,),x(tj))  is  only  used  in  the 
update  of  P(U)  but  not  in  the  gain,  since  it  was  em¬ 
pirically  shown  that  this  procedure  leads  to  an  unbiased 
estimate  of  the  state. 

4.3.  Estimation  of  u  and  Tgo 

Since  the  target  angular  velocity  term  is  embedded 
in  the  states,  u>  should  be  reconstructed  using  the  esti¬ 
mated  states.  A  simple  way  to  determine  the  value  of 

ax  ax 

u  is  to  divide  the  states  asw  =  —  or  — .  However, 

ax  alx 

since  the  expanded  state  space  is  originally  an  approxi¬ 
mated  state  space,  this  might  lead  to  numerical  errors, 
especially  when  the  higher  approximated  terms  are  used. 
By  relying  on  the  definition  of  the  vector  relation  be¬ 
tween  and  velocity  and  acceleration,  the  target  angular 
velocity  can  be  obtained  without  using  the  augmented 
states  for  approximation.  From  the  assumed  dynamics 
the  target  angular  velocity  during  its  evasive  maneuver 
is 


WT 


Vx  x  ax 

|VVI2 


(31) 


Thus,  by  using  the  state  estimates,  ui  is  constructed  as 


w  =  sign(VrzaTt  -  Vr  hr* . 


I  — I 
VT 


(32) 
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To  be  used  later  in  the  controller,  an  estimate  of  time- 
to-go,  TJ0,  is  required,  and  approximated  here  as 


Tgo  ~  ~T~  — 


\R\  \X-V\/R 


where  R  and  R  are  the  estimates  of  relative  range  and 
range  rate,  respectively,  and  the  vectors  X  and  V  are  the 
estimates  of  relative  position  and  velocity,  respectively. 


x  =  [x,y,t i,v,aTx,uTy}T 

W  =  Kt, ,  0a/JT  (36) 

Then,  the  stochastic  control  problem  is  to  find  u  which 
minimizes 


5.  Linear  Quadratic  Guidance  Law 

Based  on  the  estimated  states  and  the  estimate  of  the 
rotation  rate  constructed  from  the  estimated  states,  a 
guidance  law  can  be  mechanized.  In  the  following,  a 
stochastic  guidance  law  is  determined  which  minimizes 
a  quadratic  performance  index  subject  to  the  stochastic 
engagement  dynamics  including  the  stochastic  circular 
target  model(8)  under  the  assumption  that  states  in¬ 
cluding  the  target  states  and  the  target  rotation  rate  are 
known  perfectly.  This  assumption  simplifies  the  deriva¬ 
tion  of  the  guidance  law  enormously,  and  for  this  homing 
problem  it  is  shown  that  the  solution  to  the  stochastic 
control  problem  with  state-dependent  noise  can  be  ob¬ 
tained  in  closed  form.  The  solution  obtained  does  not 
produce  a  certainty  equivalence  controller  since  the  guid¬ 
ance  law  explicitely  depends  upon  the  system  statistics. 

Note  that  since  the  noise  in  each  cartesian  direction  is 
correlated  in  the  stochastic  circular  target  model(8),  and 
that  with  aTt  and  aTy  dynamically  coupled  through  u> 
term,  the  guidance  commands  in  the  x  and  y  direction 
cannot  be  achieved  independently.  Thus,  the  optimal 
stochastic  controller  for  circular  target  model  is  based 
on  the  minimization  of  the  performance  index 


subject  to  the  stochastic  differential  equation  with  state 
dependent  noise 

dx  =  [Ax  -f  Bu]dt  +  D(x)dd  (38) 

where 
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1  0  0  0  0  0 
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0  0  0  0 
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where  c  >  0.  Note  that 
6 

D(x)  =  £>;£>,-,  Dj€RM  (40) 

;=i 


J  —  +  y)]  +  2  J  +aAt,]^0  (34) 

subject  to  the  following  stochastic  system  of  linear  dy¬ 
namic  equations 

dx  —  udt 

dy  =  vdt 

du  =  (aTr  -aMi)dt  (35) 

de 
de 

where  Q  is  a  Brownian  motion  defined  earlier  and  £[•] 
stands  for  an  expectation  operator.  In  the  construction 
of  the  filter,  the  inherent  nonlinearity  of  the  target  model 
was  removed  by  an  expansion  of  state  variables.  How¬ 
ever,  for  the  guidance  law  formulation,  the  rotation  rate, 
w,  is  assumed  known,  although  it  must  be  constructed 
on-line  from  the  state  estimator(32). 

For  brevity  of  notation,  define  the  state  and  control 
vectors  as  follows  : 


dv  =  (aj  -  aMJdt 

darT  =  (-  JaT,  ~  vaTy)di  -  ajv 

d&Ty  =  (—  7°T,  +  waTI)dt  +  07, 


where  Xj  is  the  jth  element  of  x  and  where 

'  0  1  [0]  [0 

0  0  0 

,  D2  ,  Dzi  D4  —  Q  ,D$=  Q  ,06=  Q 

0  0-1 

0  [  1  [0 

"  _  (41) 

To  obtain  an  optimal  control  for  this  class  of  prob¬ 
lem,  dynamic  programming[13,14]  is  employed  where  the 
Hamilton-Jacobi-Bellman  equation  becomes 

0  =  J?(x,t)+Min{J°x(Ax+Bu)+hxT  A(JZx,t)x+uT  Ru}} 

(42) 

where  J°  is  the  optimal  return  function  and  the  sub¬ 
scripts  denote  partial  derivatives.  The  elements  of  the 
matrix  A  for  any  symmetric  matrix  W  is  defined  as 

Ao(^,0=<r[A(0T^;(0]  (43) 
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The  minimization  operator  in  (25)  produces 

ti  =  -R~lBTJ°x  (44) 

By  substituting  (44)  into  (42),  the  dynamic  program¬ 
ming  equation  becomes 

0  =  J°  +  J°xAx  -  l-J°xBR~lBJ°x  +  ^xTA(J^x,Q,i)x 

(45) 

The  optimization  problem  is  solved  by  explicitely 
showing  that  the  equation  above  has  a  solution  .  Asume 
J°(x,t)  =  ^ xTS(t)x ,  then 

J?  =  ^xTSx,  J°x  =  xTS ,  J°XX  =  S  (46) 

With  this  assumption,  the  dynamic  programming  equa¬ 
tion  is  satisfied  for  all  x  &  Rn  if 

S+SA+AT S+A-SBR'1  BT S  =  0,  S(tj)  =  Sj  (47) 

The  desired  optimal  controller  becomes 

u  =  -R~1BTSx  (48) 

where  S  is  the  solution  of  the  Riccati  equation  and  the 
A {S,t)ij  =  tr[DfSDj]  leads  to 

zero 

elements 

S66  -S56 
~Ss6  Sss 

The  fact  that  A  has  only  nonzero  elements  for  its 
lower-right  2x2  matrix  allows  a  tractable  closed-form 
solution.  To  see  the  characteristics  of  the  solution  in  a 
simple  manner,  matrices  are  partitioned  such  that  their 
lower-right  block  partitioned  is  a  2  x  2  matrix.  Then 
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An  A  u 

0  .4  22 
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1 A  = 


0  0 
0  S 
(50) 

where  5  is  the  nonzero  element  block  in  (49).  This  leads 
to  the  controller 


«A/X 

1 

Ri5n 

c 

B2S12 

(51) 


where  the  block  matrices  satisfies  the  decomposed  Ric¬ 
cati  equation 
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Since  the  S22  block  does  not  affect  the  block  matrices 
Sn  and  S12,  the  optimal  control  law  is  not  dependent 
on  S22 •  Therefore,  the  closed  form  optimal  guidance 
law  for  this  special  class  of  problem  can  be  obtained  by 


integrating  the  Riccati  equation  backwards  without  re¬ 
quiring  the  explicit  evaluation  of  the  A  term.  In  partic¬ 
ular,  the  stochastic  optimal  control  problem  essentially 
degenerates  to  a  deterministic  optimal  control  problem 
although  the  Ito  terms  are  retained.  The  solution  pro¬ 
cess  for  this  deterministic  control  problem,  explained  in 
detail  in  the  Appendix  A,  produces  a  guidance  law  in 
closed  form.  Note  that  the  deterministic  coefficient  A22 
includes  the  statistic  0.  Therefore,  the  resulting  con¬ 
troller  is  not  a  certainty  equivalence  controller.  The  new 
controller  becomes 
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(53) 

where  the  gains  C\  to  C4  are  an  explicit  function  of  Tg0} 
w,  and  0  as 


Ci  = 


c+  ■ 


Cl  - 


c+  : 


*  =  e*  +  (cojw'7,0  - 

c  =  c*  {-wTjoc?t-  +  (e?T-  -  cosuT,,)-  r  ~  ainuT,}^,  -“'J) 


'(f  +  w’l 


•▼here  c  = 


e?^(£+Zf)(si+u,) 


Fig.  1  is  a  block  diagram  for  an  adaptive  guidance 
scheme  for  a  homing  missile.  Note  that  guidance  gains 
are  functions  of  Tg0:  estimated  time  to  go,  the  statistic 
0  and  the  estimated  maneuver  rate  w.  Therefore,  for 
the  bearing-only  measurement  system  although  the  re¬ 
sulting  stochastic  guidance  law  is  sub-optimal  since  the 
measurements  are  nonlinear  functions  of  the  states,  the 
explicit  dependence  on  the  estimate  of  the  target  ma¬ 
neuver  rate  is  a  new  feature  which  should  help  reduce 
terminal  miss  distance. 


6.  Numerical  Simulation 

For  a  particular  engagement  scenario,  the  perfor¬ 
mance  of  the  estimator  using  the  new  target  models  and 
that  of  the  guidance  law  are  evaluated. 

6.1.  Missile  and  Target  model 

Both  target  and  missile  are  treated  as  point  masses 
and  are  considered  in  two-dimensional  reference  frames 
as  shown  in  Fig.  2.  The  missile  represents  a  highly 
maneuverable,  short  range  air-to-air  missile  with  a  max¬ 
imum  normal  acceleration  of  lOOp's.  It  is  launched  with 
a  velocity  M  =  0.9  at  a  10,000/f  altitude  with  zero 
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normal  acceleration.  After  a  0.4  sec  delay  to  clear  the 
launch  rail,  it  flies  by  the  guidance  command  provided 
by  the  linear  quadratic  guidance  law  of  Section  5  Also, 
to  compensate  for  the  aerodynamic  drag  and  propul¬ 
sion,  the  missile  is  modeled  to  have  a  known  longitu¬ 
dinal  acceleration  profile  :  ast  =  25 g's  for  t  <  2.6 sec, 
am  =  —15 g's  for  i  >  2.6sec.  The  target  model  flics 
at  a  constant  speed  of  M  =  0.9,  and  at  an  altitude  of 
10,000/t.  It  accelerates  at  9  g's  either  at  the  beginning 
or  in  the  middle  of  the  engagement.  Thus,  the  rotation 
rate  of  the  target  is  0.3  during  its  maneuver. 

Two  engagements,  considered  in  the  following  section, 
are  shown  in  Fig.  3.  With  P,  and  Rm  denoting  initial 
range  and  maneuver  onset  range,  respectively,  engage¬ 
ment  1  is  the  situation  where  the  target  maneuver  starts 
at  the  beginning,  and  for  engagement  2  the  maneuver 
starts  in  the  middle. 

C  2  Filter  Parameters  and  Initial  conditions 

Integration  of  actual  trajectories  is  performed  by  a 
fourth-order  Runge-Kutta  integrator  with  step  size  0.02 
seconds.  The  variance  for  the  angle  measurement  is  cho¬ 
sen,  as  given  in  [4],  to  be 

0  25 

Vi=a  V,  ,  Vo  =  (-^=-  +  5.625  *  10"7)/At  rad 2 

ri 

(55) 

where  R  is  range,  At  is  filter  sample  time,  and  a  is  pa¬ 
rameter  which  is  useo  in  the  simulation  indicating  dif¬ 
ferent  levels  of  sensor  accuracy. 

As  mentioned  earlier,  the  variance  for  the  pseudo¬ 
measurement  can  be  interpreted  to  show  how  strictly 
the  orthogonality  assumption  between  the  target  veloc¬ 
ity  and  acceleration  is  to  be  kept.  By  allowing  some 
acceleration  in  the  longitudinal  direction,  a  reasonable 
estimate  of  the  variance  to  be  used  can  be  given.  Sup¬ 
pose  that  the  acceleration  component  in  the  velocity  di¬ 
rection  has  a  normal  distribution  with  zero  mean.  Then 
with  probability  0.95,  a  1  g  acceleration  while  flying  with 
Vj  —  970 ft/sec  leads  to  2  a  =  3.12  *  104  [ft7 /sec3], 
where  er  is  the  standard  deviation,  which  results  in  a 
variance  V2  =  2.44  *  108  [ft7 /sec3]7. 

'Unless  otherwise  stated  the  filter  is  initialized  at 
launch  with  the  true  relative  position  and  relative  veloc¬ 
ity  component  values  assumed  obtained  from  the  launch 
aircraft.  Hence,  the  initial  values  for  the  diagonal  ele¬ 
ments  of  the  covariance  matrix  associated  with  position 
and  velocity,  Pu,  P2 2,  P33,  and  P44  are  set  to  ten  and  it 
ensures  positive  definiteness.  On  the  other  hand,  little 
knowledge  about  target  acceleration  is  assumed  to  be 
provided  at  the  beginning.  Therefore,  the  initial  values 
for  the  target  acceleration  and  expanded  states  are  zero. 
Initial  values  of  the  covariance  matrix  associated  with 
target  acceleration  is  calculated  by  resorting  to  the  def¬ 
inition  of  the  target  acceleration  at  t  =  0  given  in  (2). 
Those  covariances  are  produced  in  the  Appendix  B.  The 
target  is  expected  to  execute  a  maximum  acceleration 
turn  in  its  evasive  motion,  and  the  missile  has  no  knowl¬ 
edge  about  the  direction  of  target  rotation.  Note  that 


9  is  a  Brownian  motion  process  beginning  at  £1(0)  =  £/, 
the  expected  angle  the  target  acceleration  vector  makes 
with  respect  to  the  xr  axis  at  the  time  of  launch,  and 
a7-m>l  is  the  expected  maximum  acceleration  of  the  tar¬ 
get.  For  the  simulation  with  a  coordinate  system  having 
one  axis  perpendicular  to  the  initial  V;  direction,  0  is 
zero.  Then,  the  possible  nonzero  elements  of  the  up¬ 
per  triangular  part  of  the  initial  covariance  matrix  are 
P65(0),  P57(0),  P59(0),  Prr(0),  P79(0),  and  P99(0).  Fur¬ 
thermore,  no  information  is  available  about  the  direc¬ 
tion  of  maneuver,  and  the  possible  maximum  rotation 
rate  can  be  either  positive  or  negative.  Thus,  the  odd 
powers  of  Q  are  taken  as  zero.  This  leaves  only  Pss(0), 
P59(0),  P77(0),  P99(0)  as  the  nonzero  elements.  How¬ 
ever,  a  value  of  ten  is  assigned  to  P 661  Pgg,  and  P10  10 
to  ensure  positive  definiteness  of  the  covariance  matrix 
at  the  initial  time. 

6.4.  Filter  Results 

The  results  in  this  section  are  the  product  of  a  Monte 
Carlo  analysis  consisting  of  50  filter  runs.  Along  with  the 
miss  distance  calculations,  the  plots  of  the  estimation 
error  and  the  w  estimates  versus  time  are  mainly  con¬ 
sidered.  The  errors  are  calculated  as  [E[ex]7  +  P^y]2]1^2 
where  E[ex]  and  E[ev]  are  the  averaged  values  of  errors 
over  50  simulation  runs.  In  evaluating  the  actual  miss 
distances,  the  filter  state  estimates  are  used  in  the  guid¬ 
ance  law.  Moreover,  the  the  tracking  errors  to  be  pre¬ 
sented  below  are  based  on  the  guidance  law  in  terms  of 
the  estimated  states,  since  the  tracking  errors  were  ob¬ 
served  to  be  quite  similar  to  the  case  where  the  actual 
states  are  used  for  the  guidance  law. 

Fig.  4  represents  the  results  for  the  engagement  1 
where  the  target  maneuver  is  initiated  at  t  =  0  and  the 
pseudo-measurement  is  not  used.  It  is  shown  that  the 
estimation  improves  with  better  angle  measurements. 

When  the  auxiliary  pseudo-measurement  is  also  im¬ 
plemented  in  the  filter,  estimation  performance  im¬ 
proves  over  the  case  when  only  an  angle  measurement 
is  used.  This  is  shown  in  Fig.  5  where  again  the  target 
starts  its  acceleration  maneuver  at  the  beginning  of  the 
engagement^  =  Rm)-  At  first,  the  filter  with  the  fic¬ 
titious  measurement  seems  to  work  a  little  worse  than 
the  filter  with  angle-only  measurement.  Then,  the  fic¬ 
titious  measurement  promptly  works  as  if  it  suppressed 
or  delayed  the  filter  divergence.  Note  that  the  effect  of 
two  values  of  pseudo-noise  variance  are  shown. 

The  role  of  the  fictious  measurement  is  more  observ¬ 
able  for  engagement  2  where  the  target  maneuver  be¬ 
gins  in  the  middle  of  the  engagement  =  4000/1). 
As  plotted  in  Fig.  6,  the  filter  equipped  with  only  the 
angle  measurement  diverges  as  soon  as  the  target  ma¬ 
neuver  occurs.  On  the  other  hand,  when  the  filter  is  aug¬ 
mented  with  the  fictitious  measurement,  it  works  very 
effectively.  The  divergence  of  position  and  velocity  is  no¬ 
ticeably  suppressed,  and  the  acceleration  estimate  tend 
to  return  to  its  actual  value  from  an  instantaneous  large 
acceleration  error.  With  the  accuracy  of  the  angle  mea- 


surement  increased,  the  target  acceleralion  estimate  af¬ 
ter  the  maneuver  onset  improves  faster  than  the  filter 
that  uses  poor  angle  measurements.  This  is  shown  in 
Fig.  7. 

Performance  of  the  current  target  models  is  also  com¬ 
pared  with  the  Gauss-Markov  target  model.  The  two 
models  assume  the  same  magnitude  of  target  accelera¬ 
tion.  Note  that  in  a  Gauss-Markov  model[3][4],  A,  the 
target  maneuver  time  constant  and  IV,  the  strength  of 
the  dynamic  driving  noise  in  the  mode!,  are  two  param¬ 
eters  but  are  varied  relative  to  one  another  and  they  arc 
essestially  tuning  parameters.  However,  the  tuning  pa¬ 
rameter  is  0  in  the  new  target  model.  Along  with  the 
kinematic  constraint  incorporated  as  a  pseudomeasure¬ 
ment,  the  modified  gain  extended  Kalman  filter  is  built 
to  estimate  the  target  states,  and  the  guidance  law[4] 
is  based  on  the  Gauss-Markov  target  model.  Figs.  8- 
9  show  a  measure  of  how  well  the  filter  estimate  the 
target  states.  It  is  noted  that  the  circular  target  model 
estimates  the  target  state  better  than  the  Gauss-Markov 
model.  It  was  also  observed  during  the  simulation  that 
the  estimation  performance  of  the  Gauss-Markov  target 
model  has  been  improved  with  the  pseudomeasurment, 
and  this  is  reflected  in  the  miss-distance  caluclations  to 
be  presented.  This  is  because  the  kinematic  constraint 
increases  the  fidelity  of  the  Gauss-Markov  target  model. 

Miss  distances  have  been  calculated  on  the  basis  of 
50  runs  of  Monte  Carlo  simulations  with  an  approxi¬ 
mate  error  ±0.02  ft  due  to  subdiscretization  near  the 
final  time.  In  Table  1,  the  actual  states  are  fed  to  the 
guidance  law  in  the  Case  I,  and  the  estimated  states  and 
maneuver  rate  estimate  are  fed  to  the  guidance  law  in 
the  Case  II  and  III.  The  estimates  are  obtained  from 
angle-only  measurements  in  the  Case  II,  and  from  both 
angle  and  pseudo-measurement  in  the  Case  III.  Miss 
distance  performance  is  tested  as  more  noise  is  intro¬ 
duced  into  the  measurement  and  then  into  the  dynam¬ 
ics.  Table  1  indicates  that  much  of  improvement  comes 
:rcm  the  circular  target  model  with  additional  improve¬ 
ments  achieved  from  the  kinematic  constraint.  In  addi 
tion,  miss  distance  has  been  improved  by  using  the  an¬ 
gle  and  pseudo-measurement,  especially  as  the  process 
noise  power  spectral  density  ©  in  the  state  dependent 
noise  term  decreases.  Calculation  of  miss  distance  with 
the  Gauss-Markov  target  model  also  indicates  that  sig¬ 
nificant  improvement  is  obtained  in  the  Guass-Markov 
model  using  the  kinematic  constraint.  For  the  particular 
scenario  chosen  here,  circular  target  model  augmented 
by  pseudo-measurement  out  performs  the  Gauss-Markov 
target  model. 

7.  Conclusions 

ii.e  orthogonality  between  the  target  acceleration  and 
velocity  vectors  is  a  typical  characteristic  of  the  target 
of  an  air-to-air  missile  and  it  is  utilized  in  the  develop¬ 


ment  of  a  new  stochastic  target  acceleration  model  for 
the  homing  missile  problem.  In  addition,  this  character- 
istic  is  also  implemented  in  the  form  of  an  augmented 
pseudo-measurement.  A  guidance  law  that  minimizes 
a  quadratic  performance  index  subject  to  the  stochas¬ 
tic  engagement  dynamics  is  determined  in  closed  form 
where  the  gains  are  an  explicit  function  of  the  estimated 
target  maneuver  rate  and  time  to  go.  Preliminary  results 
for  the  two-dimensional  case  indicates  that  the  circular 
target  model  is  able  to  produce  a  reliable  estimate  in  the 
homing  missile  engagement.  When  it  is  augmented  bv 
the  fictitious  measurement,  tl  e  modified  gain  extended 
Kalman  filter  using  the  proposed  target  model  results 
in  the  significant  enhancement  of  target  state  estima¬ 
tion.  The  kinematic  constraint  also  leads  to  the  signifi¬ 
cant  improvement  in  miss  distance  performance  for  the 
Gauss-Markov  target  model.  Comparisons  of  the  cur¬ 
rent  target  models  over  the  Gauss-Markov  target  model 
show  that  a  significant  improvement  is  gained  in  target 
state  estimation  and  miss  distance. 
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Appendix  A 

Linear  quadratic  guidance  law  for 
deterministic  circular  target  model 
In  the  following,  the  optimal  deterministic  guidance 
law  for  linear  quadratic  problem  is  sought  for  the  current 
circular  target  model  filter.  The  det  mmi  nstic  optimal 
solution  can  be  obtained  by  solving  the  Riecati  equation 
without  tne  A  term  via  transition  matrix  approach,  but 
the  use  of  Euler-Lagrangc  equation  seems  simpler  for 
this  case. 

The  problem  is  to  minimize  the  performance  index 

J  =  ^ +  vj]  +2 

subject  to  the  following  linear  dynamic  system 

x  =  u 
y  =  v 


where  A  ,i  —  1 , ..., 6  is  a  Lagrange  multiplier.  The 
Euler-Lagrange  equations  for  A,-  are 

Ai  =  9,  A2  =  0,  —A3  =  A!,  -A«  =  A2 

w'herc-  the  optimal  control  satisfies  the  optimality  condi- 


Hnally,  the  Euler-Lagrange  equations  with  the  natural 
boundary  conditions  yield 

Ai  =  x j ,  A2  =  y/ ,  A3  —  xjTgot  A4  —  yjTgo 

which  gives  the  control 

&M,  =  x(tf)T}0/c,  =  y{tj)Tgo/c 

where  T}J  is  the  time-togo  of  missile  to  intercept  the 
target  and  c  is  the  guidance  law  design  parameter.  In 
order  to  get  the  guidance  law  in  terms  of  the  current 
states,  the  underlining  dynamics  is  integrated  backward 
from  tj  to  t.  Succesive  integrations  of  state  differential 
equations  yield 

or,  =  ecja>7’,„e?r»-aT,(!/)  +  ainw7',.e?T'-a7>(*/) 
arw  =  -■«nw7',.e^r'*or.(t/)  +  co5-;r,ae?r>*oT>(t;) 

u  =  ^r.V((/)  +  “(</) 

1  © 

+1J— — jIttO  ~  c°suTlo'7T~)  -wjinw7',„c?T-]aT.(l/) 

+  1T^7!“u’(1  ~  eMu,r»»e?r")  -  f  ■ «nu>7‘,.e?Ji*)la7-,(j/) 

"  =  +  «( </) 

“  C03uT,^T”)  +  |iiW7',ce?T>*]oT.((/) 

+  «!— -cotuT„t$T~)  -  u>iinaiT„ti7-)]aT,(ti) 


u  =  ax,  -  a 
V  =  ax,  ~  aAf, 

0 

aX,  =  -  yar,  “waT, 

0 

ox,  =  ~  jOTi  +  war. 

This  linear  system  of  dynamics  stems  from  taking  Ito 
derivative  of  the  corresponding  nonlinear  stochastic  tar¬ 
get  model  (6).  The  w  is  the  angular  rate  of  target  maneu¬ 
ver  which  is  handled  as  a  known  constant  in  the  deriva¬ 
tions.  In  the  actual  mechanization  of  guidance  command 
the  value  of  Q  constructed  from  the  estimated  states  are 
used. 

The  variational  Hamiltonian  and  the  augmented  end¬ 
point  function  are  given  by 

H  =  ca\f*  -(-  ca.\t*  +  Aiu  +  A2v  +  A3(a7r  —  a^f,) 

+  ^a(ox,  ~  +  A5(-—  -  war,)  +  A6(-—  +waxj 

G  =  \ +yj) 


x  =  d--^ )*(«;) -r,„u(r;) 

+ 7^~Tj3inuTt>e^arAtf) 

+  E^flu,:r«» +  (Srfbj(1  “  co"7>?7V*) 

+  s’- - p\~wTio  ~  7Sr~xrO  -  cosu>T„crT“) 

TT“'  (t+“J 

(—  —  u>2) 

w0  *T  , 

+  nn— — Jinwr,.eTr-|or,(</) 
lT+u) 


The  final  states  being  expressed  in  terms  of  the  current 
states  via  6x6  matrix  inversion, the  optimal  guidance  law 
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is  obtained  as  equation  (51).  As  expected  from  dynamic 
coupling  in  the  target  acceleration  model,  guidance  com¬ 
mands  in  each  channel  are  the  function  of  acceleration 
components  in  both  x  and  y  directions. 

Appendix  B 

State  and  Error  variance  associated  with  target  acceleration 


Since  the  initial  values  for  the  state  estimates  associ¬ 
ated  with  target  acceleration  are  set  to  zero,  the  state 
and  error  variances  are  computed  with  the  aid  of  ex¬ 
pected  values  of  trigonometric  functions  such  as  Table  1.  Comparison  of  miss  distance 


/•  +  OO 

£[cos2  6}  =  /  cos ^  0[{6)d6  with  p(6) 


By  using  standard  manipulation,  the  expected  value  of 
the  cos'  6  is 


f?[cos20]  =  ^[1  +  cos26e~‘e'], 

and  in  the  same  manner 

£[sin2  0]  =  ^[1  H -  cos20e-2e‘], 

£[cosf?sin0]  =  ^sin20e_2e‘. 

This  veilds  the  initial  conditions  for  the  state  and  error 
variances  associated  with  target  acceleration  as  follows. 


/VO)  =  -Y«(0)  =  «iL.  *  A 
/V(0)  =  Xst(O)  =  *  /', 

/W 0)  =  A'„(0)  =  ■  / 1, 

/WO)  =  Af«f0)  =  •/ 3, 

P n(0)  —  -AW0)  —  ar  , , **•  *  / 3, 
/eio(0)  =  Agio(O)  =  •  /3 1 

Pr,(0)  =  X„(0)  =  .  / 2, 

Pno(O)  =  Afr.ofO)  =  .  / 2, 

/W0)  =  A'w(0)  =  4„„w3*/2, 
/WO)  =  X»(0)  =  *  / 1, 

P ioio(0)  =  Aicio(O)  —  *  /3 


/’WO)-  XWO)  =  «!-.„*  /2, 

/WO)  =  Au(0)  =  .  /2, 

/iio(0)  =  A'jio(O)  -  *  / 2. 

P67(0)  =  X,,(0)  =  4.„**/2. 

/WO)  =  A«s(0)  =  4,„W  .  /2, 
Pn(0)  —  AVr(O)  —  *  /l, 

/WO)  =  A„(0)  =  »  /l, 

P«(0)  =  A„(0)  =  a\mJ*  •  /3, 

/»io(0)  —  Ajio(O)  —  *  /3 > 

/sio(0)  —  AsiofO)  ==  or*,,,**’4  * 


(1  -  a.  ’5) 


where  /I  = 


(1  +  eos20) 


/  2  = 


sin20 


/3  = 


Statistic 

Case 

Range ,  (ft) 
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Figure  1:  Block  diagram  of  homing  missile  guidance 
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Figure  4:  Circular  target  model  with  different 
Vl’s(Engagement  1) 


Figure  2:  Inertial  reference  frame  for  missile  and  target 
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Figure  5:  Circular  target  model  with  pseudo- 

measurement(Engagement  I) 


Figure  3:  Typical  missile  target  trajectories 
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Figure  6:  Circular  target  model  with  pseudo- 

measurement(Engagement  2) 
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Figure  7:  Circular  target  model  with  pseudo- 

measurement(  Engagement  2) 
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Figure  8:  Circular  target  model  with  pseudo- 

measurement(Engagement  1) 


TIME/SEC)  TIME<SEC) 


Figure  9:  Circular  target  model  with  pseudo- 

measurement(Engagement  2) 
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